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Abstract 



This paper shows that the ordinary Brillouin spectrum peaks asso- 
ciated to scattered radiation off acoustic modes in a fluid suffer a shift 
^ ' in their values due to gravitational effects. The approach is based in 

■ the ordinary linearized Navier-Stokes equations for a fluid coupled to 

a Newtonian gravitational potential. The formalism leads to a disper- 
sion relation that contains both gravitational and dissipative effects. 
It is also shown that the Brillouin peaks tend to condense into a single 
peak when the fluid modes approach the critical Jeans wave number. 



1 Introduction 

Dynamic light scattering theory PQ-|2] shows that photons interact with the 
acoustic modes in a fluid. This scattering arises from sound waves associ- 
ated to statistical fluctuations in the fluid's density and temperature. Early 



1 



theoretical work in this subject, as well as the measurements confirming the 
existence of the Brillouin-Rayleigh (BR) spectra are given in Refs. [l]-[3]. 
The first theoretical approach to the subject is a lapidal note published by 
Landau and Placzek back in 1936 [4] which has been the inspiration of most 
of the work done on this problem. Recently, these ideas have been extended 
into the realm of astrophysical and cosmological contexts in relation with the 
discussion of the old problems on gravitational instabilities and Jeans' mass 
number [5] jH] • The recognized fact that Cosmic Microwave Background Ra- 
diation (CMBR) photons interact with acoustic fluid modes in certain stages 
of the universe evolution suggests that Brillouin scattering is potentially im- 
portant to CMBR physics. 

On the other hand, entropy production plays an important role in re- 
alistic descriptions of thermodynamical processes, and some discussion of 
its physical sources may yield interesting results in the analysis of the time 
evolution of statistical fluctuations in astrophysical systems [TT]. Recent 
work JU] shows indeed that the critical Jeans wave number for gravita- 
tional collapse is slightly modified by dissipative effects such as viscosity. 
The formalism used in that work is also useful to obtain structure factors, 
namely scattering laws that describe line broadenings and frequency shifts 
associated to the interaction between the incoming radiation and the fluid 
modes. The mathematical tool in this case is the so-called density-density 
correlation function, taken in the frequency (Fourier) domain which is ob- 
tained by standard techniques [2|-[3]. The influence of the gravitational field 
in the Brillouin spectrum will become apparent from the analytic form of the 
structure factor. 

In this paper we wish to discuss two results. Firstly, that the presence of a 
gravitational field produces a shift in the effective frequencies of the Brillouin 
doublet in a simple fluid. The shift could be significant for large enough 
densities such as those of matter in early stages of the Universe. Secondly, the 
Brillouin peaks tend to condense into a single peak as the critical Jeans wave 
number is approached. This implies that at certain densities, the thermal 
fluctuations and the "mechanical dissipative modes" behave more or less in 
the same way, a fact which, if detectable, should become manifest in the 
anisotropies of the observed CMBR. The significance of this fact has not yet 
been explored. 

This paper is thus divided as follows. Section two is dedicated to a review 
of the linearized form the Navier-Stokes equations in the presence of a grav- 
itational field. Section three is devoted to the derivation of the scattering 
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law and its comparison with the ordinary BR spectrum. A brief discussion 
of the implications of the results here obtained is included in section four. 



2 Linearized Navier-Stokes system in the pres- 
ence of a gravitational field 

The standard system describing fluctuations of the thermohydrodynamical 
variables in a simple non-reacting fluid consisting of particles of mass m in the 
presence of a gravitational field, is derived by following the tenets of Linear 
Irreversible Thermodynamics (LIT) [2] JUj. Indeed,the so-called Navier- 
Stokes-Fourier equations of hydrodynamics for a simple fluid arise from the 
structure of what is now called LIT by supplementing the two conservation 
equations for mass and momentum, respectively and the balance equation 
for the internal energy, with additional information. Indeed, these equations 
read 0-13 ■ 

D e dQi _ pi7 . dm 
' Dt dxi " dxi 1 } 

Here, p(a?,t), u l {x^,i) and e(x J ',t) are the local density, velocity and in- 
ternal energy, respectively, S 4 - 7 the momentum current (or stress tensor) and 
Qi the heat flux. All indices run from 1 to 3. In general, for isotropic flu- 
ids, H lJ = p5 l i + r* 3 ', where p is the local hydrostatic pressure and r*- 7 the 
viscous tensor (S 1 ^ is the unitary dyadic). Notice that Eqs. (PE1) contain 
fifteen unknowns (including the pressure) and there are only five equations, 
so the system is not well determined. If we arbitrarily choose to describe 
the states of the fluid through the set of variables p(x^,t), u l (x^,t), T(x^,t), 
where T is the local temperature, we need nine dynamic equations of state 
(or constitutive equations) relating and Qi to the state variables plus 
two local equations of state p = p(p,T) and e = e(p,T). According to the 
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tenets of LIT we choose the so-called linear constitutive laws, namely 



T 



-2na l] -C05 13 



(4) 




(5) 



with 77 and £, the shear and bulk viscosities, respectively, a*- 7 the symmetrical 
traceless part of the velocity gradient, = jfi and k being the thermal 
conductivity. 

Eqs. (@J) and (0) are the well known constitutive equations of Navier- 
Stokes and Fourier, respectively. Substitution of these equations into Eqs. 
(PCH) yields a set of second order in space, first order in time non-linear 
coupled set of partial differential equations for the chosen variables p, u % and 
T. This set, which the reader may seek in the literature [21,0-111, is the 
so called Navier-Stokes-Fourier system of hydrodynamic equations. The non- 
linearities appearing in such equations have two sources, the inertial terms 
u 1 ^; arising from the hydrodynamic time derivatives, plus quadratic terms 
in the gradients of velocity arising from Eqs. (j3J) and ©. Moreover, it should 
be mentioned that this set of equations is consistent with the second law of 
thermodynamics, the Clausius uncompensated heat , or entropy production, 
is strictly positive definite. 

Nevertheless, for the purpose of this paper, this set of equations is too 
complicated. In order to deal with fluctuations around the equilibrium state, 
one assures that for any of two state variables, call them X(x^,t) one can 
write that, 



where X a is the equilibrium value of X and SX the corresponding fluctuation. 
Neglecting all terms of order (5X(x\ t)) 2 and higher in the NSF non-linear 
set one finds the linearized NSF equations of hydrodynamics 0-[Zl, 



X(x j ,t) = X + SX(x j ,t) 



(6) 



d_ 

dt 



(5p)+ P 6 = 



(7) 



du k 



1 d 



(Sp) 



(3 d 
kt dx k 



(5T) + 2 V V 2 u k -^ V -,)^:(e) + F k (8) 



Po 



dt 



p k t dx k 



d_ 

ot 



(6T) = D T V 2 5T 




e 



(9) 
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since u ko = 0,u k = 5u k and ^ = (§g) Po , k t = p (§^)r and D T = is 
the thermal diffusivity. C p and C v are the specific heats at constant pressure 
and constant volume, respectively. 

Note now that from Eq. (jHJ), u kTrans = Rot (u k ) uncouples from the hy- 
drodynamic modes, whence, one formally arrives at the desired set, namely: 

IL( 6p) + Po e = o (io) 

Po ^- = - -^—W 2 6p - ^-V 2 5T + D V V 2 56 - V 2 ^ (11) 
ot p Kt Kt 

I- (ST) = D T V 2 5T - -^-9 (12) 

assuming that F k = — J^, where p is the gravitational potential. Also, 
D v = (|?7 + Now, we notice that (5p) = —p o so we can reduce 

this set to only two equations, Eq.ffTTj) and Eq.(12). We now perform a few 
minor transformations introducing the speed of sound of the fluid C 2 Q = 
through the relationship 

K T = 7K S = (13) 

where k s is the adiabatic compressibility (j^j whence, = • This leaves 
us finally with the set: 

- ^T- + - V2 ^ + — V 2 <5T + D v V 2 d (Sp)) + Po V 2 5 V = (14) 
ot 7 7 Ot 

| (iT )_ Dr V^-l^|(V) = (15) 

where we have used that C v — C v = ^r- 2 -- 

Eqs. (11411 form a set of coupled equations for the density and tempera- 
ture fluctuations in the fluid under the action of a conservative force whose 
nature need not to be specified for the time being. They are the basis for 
studying the properties of the time correlation functions of thermodynamic 
fluctuations. Those of the density will be of particular interest here. 
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3 Solution to the hydrodynamic equations 



The results derived in the previous section are far from being new. Aside 
from the term V 2 which arises from the presence of an external conservative 
force, they are identical to the ones that have been widely discussed in the 
literature. The question here is if the gravitational potential introduces any 
substantial modification in the correlation functions for the thermodynamic 
fluctuations. To examine this possibility we recall that, if we are consid- 
ering only the fluctuations, the gravitational potential satisfies the Poisson 
equation: 

V 2 5<p = -4nG5p (16) 

Now, equation (JT5|l reads 

- + ^V 2 5p + ^V 2 5T + AV 2 (| (dp)) - AnG Po 5p = (17) 

The introduction of the Poisson equation links fluctuations in the gravita- 
tional potential with density fluctuations. The solution to Eqs. (fT5j) and (JT7|) 
proceeds in the standard fashion. We reduce them to a set of algebraic 
equations by taking their Laplace-Fourier transform, choose to set the static 
temperature fluctuations equal to zero, and eliminate the temperature lead- 
ing to an equation for 5p (k, s^j , which is the ratio of two polynomials in s. 
In fact, one gets that 

8p(k, s) _ s 2 + (A + Dt)k 2 s + k*D v D t + C k 2 (l - i) 

6p(k, 0) ~ (s + A fc 2 )(s 2 + D v k?s + C 2 P - AttG Po ) (18) 

To compute 5p (k,tj one must take the inverse Laplace transform of the 
former quantity, which demands the knowledge of the roots of the denomi- 
nator, which is a cubic equation in s (dispersion equation). 

The analysis leading to Eq. (jl8|) clearly points out the fact that the solu- 
tion to the cubic equation giving rise to the poles of the function 5p(k,s), 
is exact. This is an improvement over the current version in the literature 
[I] asserting that the roots are only approximate to order k 2 . There 
are two immediate implications of this result, namely, Rayleigh's peak is not 
affected by the gravitational field within the Navier-Stokes regime and the 
Jeans number, thoroughly discussed in Ref. [TUj cannot be affected by thermal 
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dissipation, the thermal conductivity will never appear in its definition, as 
will be pointed out below. The former result appears to be in contradiction 
with previous work stating that Rayleigh's peak is modified by a constant 
gravitational force jH]. Yet, notice that this is a k A effect, which is beyond 
the Navier-Stokes domain. Work along this line with the Burnett and Su- 
perBurnett equations is in progress. So, we shall leave the discussion of this 
subtlety for the future. 

Returning to Eq. (fTHj). the two roots of the quadratic equation in the 
denominator are: 



S l,2 



D v k 2 



±i 



Clk 2 



Dlk 4 



nl/2 



(19) 



If 4irGp = and viscosity dominates over the term in k 2 , this result reduces 
to the one giving rise to the standard Brillouin peaks which correspond to 
density fluctuations of the type |2j [Zj, 



5p (k, t) = 5p (k, 0) -e- D " kH Coa [C a kt] 



(20) 



which are the acoustic modes damped by the Stokes-Kirchhoff factor D v . 

On the other hand, if 4nGp ^ 0, the threshold value for k distinguishing 
between damped oscillations and growing modes is given by 



(C 2 k 2 - AixGpo 



D 2 k A 







(21) 



or, 



2CI 
D 2 



1 ± 



\ 



^Gpo£l 



(22) 



Eq. (j22j) is a generalization of Jeans wave number when dissipative effects 
due to viscosity are non- negligible, and is the main result of the paper. We 
can note that, if 47r< ^f 4 ° D " -C 1, and taking the — sign for the square root, we 
have: 



or 



k 



AnGp 



C 2 



(23) 
(24) 
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This is the square value of the Jeans wave number. It was derived by Jeans in 
1902 and rederived in many other waves by several authors. Here we simply 
show that it is almost a trivial consequence of LIT. Jeans wave number has 
been used by cosmologists to estimate the minimum mass required to form 
a galaxy. This question has been discussed at length in Ref. 

Coming back to our main result, Eq. (|18|) . it is clear that damped acoustic 
waves are associated to the case in which the dispersion equation has complex 
roots. The corresponding expression for the density fluctuations in the (k, t) 
space, neglecting the last term in the square root of Eq. ()19|). is then given 
by: 



(C 2 k 2 -AnG Po ) 1 ' 2 t 



6p(k, t) = 5p(k, 0) 1 1 1 - -J e~ UT k 1 + -e~ v k 'Cos 

1 (25) 
Eq. (f23j) represents damped acoustic waves propagation with an effective 

1 /2 

frequency (C 2 k 2 — 4nGp ) . When gravity is negligible, one recovers the 
usual frequency u> = ±C Q k. Notice, however, that although G is a small 
number, it is multiplied by the density p . In fact, it may turn that for a 
critical value of p Q the damped modes are replaced by growing modes (Jeans 
instability). 



4 Structure factor and modified Brillouin spec- 
trum 

Following the standard dynamic radiation scattering theory [2] -[3], one may 
then construct the density- density correlation function from Eq. (j25j) . This 
leads to the expression: 



(6p*(k,0)6p(k,t)) 
(<5p*(k,0)5p(k,0)> 



e-° Tkt + 



- Vk2t Cos 
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(C 2 k 2 -ATtG Po ) l,2 t 



(26) 

where the brackets < > indicate an average over an equilibrium ensemble. 

The Brillouin specific intensity of the scattered light due to its interaction 
with the acoustic modes of a fluid is obtained by Fourier's transform of Eq. 
The ensuing expression reads: 
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5 /¥ ,(k, w ) = J[(l-i) ;3 ^ 5I + i 



Tk 2 



7 ' lu 2 —D^ fc 4 7 V ( w _(c2fc2-47rGp ) 1/2 ) +r2 fe4 

rfc 2 

(w+(C 2 fc2_ 47rGpo )l/2-)2 +r2A;4 , 



r,2 ; » (27) 



The last two terms in Eq. (}2*Tj) correspond to the Brillouin doublet and 
deserve some attention. This doublet reflects a shift in frequency governed 
by the speed of sound of the fluid, the wave number of the incoming radia- 
tion, and the effect of gravity. This last effect can be roughly estimated by 
performing a binomial expansion of the square root included in the effective 
wave frequency. Indeed, let u p , the modified frequency be defined as: 

u P = ± (C 2 k 2 - 4ttG Po ) 1/2 ~ C a k (l T (28) 
to first order in Gp Q . Thus, if we now call |Ai/ p | = \uj — u p \ we see that 

\*>,\ = ^ (29) 
and the relative change in the frequency Brillouin peaks is, roughly: 

\Av p \ ( Gp c 2 \ 1 



C a k \Ait 2 C 2 



(30) 



if using k — -y — ^-za Written in this way the frequency exhibits a much sub- 
tler dependency on other variables such as the temperature, which implicitly 
appears in C Q . 



5 Final remarks 

A rather interesting result appears in Eq. JZZD if k = kj = ^3^, the Jeans 
wave number. In that case, 5p(k, t) is a simple decaying exponential with 
time, and the Brillouin-Rayleigh spectra reduces to the ordinary Rayleigh 
peak superposed to a similar one arising from the collapse of the Brillouin 
peaks, suggesting that the Stokes-Kirchhoff dissipative factor behaves in a 
way similar to the thermal fluctuations. This collapse of the Brillouin peaks 
would cause a line broadening in scattering situations and would reflect a 
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competition between entropy fluctuations (D?) and mechanical fluctuations 
(D v ). Whether or not this result is of any relevance at all remains to be tested. 

A final remark concerns Eq. (|29|) . It is well-known that for typical electro- 
magnetic radiation, the frequencies are rather high, even radio frequencies. 
This is not the case for gravitational waves, whose wavelengths are very long. 
One could then speculate about the possibility of having gravitational radi- 
ation Brillouin scattered by cosmological matter. This would enhance the 
relative change in frequency of the Brillouin peaks. 

This work has been supported by CONACyT (Mexico), project 41081-F. 
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